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ABSTRACT 

We report observations of quasi-periodic pulsations (QPPs) during the X2.2 flare of 15 Febru- 
ary 2011, observed simultaneously in several wavebands. We focus on fluctuations on time scale 1-30 s 
and find different time lags between different wavebands. During the impulsive phase, the Reuven 
Ramaty High Energy Solar Spectroscopic Imager (RHESSI) channels in the range 25-100 keV lead 
all the other channels. They are followed by the Nobeyama RadioPolarimeters at 9 and 17 GHz and 
the Extreme Ultra- Violet (EUV) channels of the Euv Spectrophotometer (ESP) onboard the Solar 
Dynamic Observatory (SDO). The Zirconium and Aluminum filter channels of the Large Yield Ra- 
diometer (LYRA) onboard the Project for On-Board Autonomy (PROBA2) satellite and the SXR 
channel of ESP follow. The largest lags occur in observations from the Geostationary Operational 
Environmental Satellite (GOES), where the channel at 1-8 A leads the 0.5-4 A channel by several 
seconds. The time lags between the first and last channels is up to ~ 9 s. We identified at least two 
distinct time intervals during the flare impulsive phase, during which the QPPs were associated with 
two different sources in the Nobeyama RadioHeliograph at 17 GHz. The radio as well as the hard 
X-ray channels showed different lags during these two intervals. To our knowledge, this is the first 
time that time lags are reported between EUV and SXR fluctuations on these time scales. We discuss 
possible emission mechanisms and interpretations, including flare electron trapping. 
Subject headings: Waves — Sun: flares — Sun: oscillations — Sun: radio radiation — Sun: UV 
radiation — Sun: X-rays, gamma rays 



1. INTRODUCTION 

Quasi-periodic pulsations (QPPs) have been observed 
during solar flares for many y ears in many wavelengths, 
from radio to hard X-rays (Nakariakov and Mclnikov 
l2009f ) and even gamma rays ((Nakariakov ct al. 2010) 7 . 
The observed periods range from fractions of a sec- 
ond to several minutes. Recent work reflects this 
divers ity in time s c ales a n d wavelengths of observa - 
tion: iFoullon et all (120101): iKuprivanova et all (|2010l ): 
iReznikova and Shibasaki (2 011[). 

iNakariakov and Melnikovl (|2009f ) divide the current 
interpretations into two classes: periodic load/unload 
mechanisms of non-thermal electrons produced during 
the flare (i.e. an intrinsic property of the reconnection 
mechanism that leads to a quasi-periodic behavior) or 
modulation of the electron beam or of parameters of 
the emi tting plasma by magnet o-hydrodynamic (MHD) 
waves. iFleishman et al.l ([2008), for example, interpret 
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their observations in the framewo rk of quasi-period ic in- 
jection of fast electrons, while for lAsai et al.l (|2001h . the 
inject ion is modulated by an oscillation of the flaring 
loop. Ingl is et al.l (|2008f) favor an MHD wave, but do not 
exclude a quasi-periodic injection of the fast electrons. 

In this Letter we investigate this question by compar- 
ing the signal observed in different wavebands during a 
solar flare. We pay particular attention to the time de- 
lays and focus on short time scales fluctuations (< 30 s). 

2. OBSERVATIONS 

Many instruments observed the X2.2 flare that started 
at 01:44 UT on 1 5 February 2011 in AR 11158 
(|Schrijver et al.ll2011h . We focus here on the instruments 
that provide high time cadence (cf Tabl e Ht . 

Th e Large Yield Rad i omete r (LYRA; IHochedcz et al.1 
120061: iBenmoussa et al.l I2009D on the Project for On- 
Board Autonomy (PROBA2) spacecraft observes at up 
to 100 Hz in four channels. We focus here on the Zirco- 
nium (Zr) and Aluminum (Al) channels, which both in- 
clude contributions from soft X-rays (SXR) and Extreme 
Ultraviolet (EUV); see Table Q] for details. LYRA's two 
UV channels, Lyman-a and Herzberg-continuum, con- 
tained no significant fluctuations, probably due to the 
degradation of the filters of the unit in use at the time 
of observation (Domi nique et al1l2012h . 

The Nobeyama RadioPolarimeters (NoRP) measured 
for this event solar fluxes at 1, 2, 3, 9, 17 and 35 GHz. 
We also used reconstructed images from the Nobeyama 
RadioHeliograph (NoRH) at 17 GHz 8 . 

The Euv Spectrophotometer (ESP; iDidkovskv et al.l 

8 available at http : //solar . nro . nao . ac . jp/norh/images/event/20110215_01 
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I2012D is part of the Extreme ul t ravio let Variabil- 
ity Experiment (EVE; iWoods et al.l I2012D instrument 
suite onboard the Solar Dynamic Observatory (SDO; 
iPesnell et al.1120 12). It comprises five channels: one SXR 
and four EUV channels. The 18 nm waveband contains 
line emission that mainly comes from a plasma around 1 - 
2 MK. The 26 nm and 30.4 nm channels mainly contain 
He II lines. We subtracted an offset to correct an artifi- 
cial jump in instrument bias that appeared beginning at 
02:00 UT in all these channels. (We did not use EVE's 
36 nm channel because of a technical problem with that 
channel.) 

The Geostationary Operational Environmental Satel- 
lite (GOES) provides SXR fluxes in a "short-wavelength" 
bandpass (0.5-4 A) and a "long-wavelength" bandpass 
(1-8 A). 

The Reuven Ra maty High Energ y Solar Spectroscopic 
Imager (RHESSI; iLin et all I2002T ) observes from 3 to 
20000 keV. We restricted our analysis up to 300 keV, 
the emission being absent above that limit. 

We also used the A tmospheric Imaging Assembly 
(AIA; ILemenet al~ll2011l) EUV imager onboard SDO for 
context identification. 

3. RESULTS 

Panel A in Figure Q] presents the light curves we ob- 
served in several different channels. The Nobeyama 
17 GHz channel, like the RHESSI channel (not included), 
shows the impulsive phase of the flare. The other chan- 
nels peak with different delays, depending on their re- 
spective waveband and temperature response. Panel B 
shows several irradiance curves with the large-scale trend 
removed, which we achieve by subtracting the corre- 
sponding signal smoothed with a 20-s boxcar. For sim- 
plicity we omit channels from each instrument that fea- 
ture similar properties. 

Clearly, short time-scale fluctuations appear during the 
impulsive phase for each of the plotted channels. Some 
fluctuations of smaller amplitude are also present during 
the declining phase, for example, in the SXR channels. 
In most cases, the fluctuations have a small amplitude 
relative to the overall increase in signal during the flare 
(Tabic P. 

Panel C shows an enlarged view where we superim- 
pose every channel (black) over the ESP SXR channel 
(red) . Although the curves of total irradiance have dif- 
ferent shapes, especially the radio channel, the short term 
fluctuations (< 20 s) are strikingly similar. 

In most cases, time delays appear between similar fea- 
tures. We investigated further by measuring the time 
delay that provides the maximum correlation between a 
channel and the ESP SXR channel, which we consider 
our reference channel; best correlation coefficients are 
shown in Table Q] We used linear interpolation to match 
the time cadence of the detrended signals (spline interpo- 
lation produced similar results) . We also smoothed each 
signal using a 1-s boxcar to reduce the effect of noise. 
This smoothing has a negligible effect, but the curves of 
correlation coefficient as a function of time lag are nois- 
ier and contain some spurious local maxima without it. 
Additionally, we verified that the best lags were not sig- 
nificantly different when using smoothing boxcars with 
widths from 5 to 30 s to remove the large-scale trend. 
We also cross- verified the optimal correlation lags by us- 



ing several different channels for reference. 

It is clear that some local extrema in the RHESSI chan- 
nels are nearly in anti-phase with those in ESP SXR, 
while some other extrema are in phase. As most extrema 
in the RHESSI channels from 3 to 100 keV coincide, the 
effect of Poisson noise can be ruled out. Nevertheless, 
it is difficult to distinguish between the effect of signal 
integration over the 4-s rotation of the spacecraft, in- 
strumental noise, and real differences with the ESP SXR 
channel due, for example, to different contributing pro- 
cesses. Such differences could also be the result of the 
fact that soft and hard X-rays (HXR) originate in differ- 
ent parts of the flaring region. 

Figure[2]synthesizes the delay times with respect to the 
ESP 0-7 nm channel for channels whose best correlation 
coefficients were larger than 0.4. There are three intervals 
over which the fluctuations exhibit clearly distinguish- 
able behavior. Interval la lasts from 01:48:00 to 01:50:15; 
interval lb lasts from 01:50:15 to 01:59:00. These inter- 
vals correspond to the early and late impulsive phase of 
the flare. Interval 2, from 01:59:00 to 02:10:00, corre- 
sponds to the portion of the declining phase where fluc- 
tuations are significant. 

The main difference between intervals la and lb is that 
some radio channels lead the reference channel during in- 
terval la, but do not during interval lb. In fact, consid- 
ering possible offsets of the various instruments' clocks 
(believed to be < 1 s), it is possible that the radio and 
ESP 0-7 nm channels are actually in phase during in- 
terval lb. RHESSI time series also contain different lags 
during these two intervals as well. Thus we note that this 
behavior occurs particularly in wavebands characteristic 
of the impulsive phase. 

All the other channels retain roughly the same lags 
throughout the impulsive phase. The three EUV chan- 
nels are nearly all in phase and lead all the SXR channels. 
Within the EUV and SXR channels we find, roughly, that 
harder spectral components contain larger lags. Addi- 
tionally, both LYRA channels slightly lead the ESP refer- 
ence channel, which is likely due to combination of EUV 
and SXR emission that contributes to these wavebands. 

During the declining phase (Interval 2), fluctuations 
appear in all SXR channels (but sometimes with correla- 
tion coefficient slightly smaller than 0.5), but not in the 
HXR, EUV or radio channels. 

This method used above disregards the frequency com- 
position of the signal and, moreover, only responds to 
short time-scale fluctations. Computing the wavelet 
coherence between the sign als, using the method of 
iTorrence and Websterl (|1999f ) 9 may return additional in- 
formation and can be performed without removing the 
long-term trend, but this method also has limitations. In 
particular, no useful information can be recovered for sig- 
nals that are very similar on long time-scales, such a flare 
pulse, because all periodicities are found to be coherent 
and therefore cannot be distinguished. Nonetheless, it 
did prove useful for some of our data. 

Figure [3] shows the wavelet coherence between the ESP 
SXR and 18 nm channels. The signals have been first 
normalized by subtracting the average and then dividing 
by the standard deviation (top panel). The detrended 
signals are only shown in the middle panel for comparison 

9 software available at http : //atoc . Colorado . edu/research/wavelets/ 
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Fig. 1. — Normalized irradiances in different instruments during the 15 February 2011 X2.2 flare (Panel A). To emphasize the short-period 
fluctuations, Panel B shows the same observations detrended by subtracting the signal smoothed using a 20-s boxcar (in units of their 
standard deviation). Panel C shows a close-up view over the interval delimited by the vertical dotted lines in Panel B. The ESP 0—7 nm 
light curve (red) is over-plotted on each curve for comparison. The tickmarks are drawn at ±2 standard deviations of the signals. The 
artifacts in the LYRA Zr curve around 01:56 and 02:20 are due to PROBA2 spacecraft maneuvers. 
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Fig. 2. — Time delays for the detrended irradiances calculated with respect to the detrended irradiance of the ESP 0-7 nm channel (a 
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Fig. 3. — Normalized light curves of the SXR and EUV 18 nm channels of SDO/EVE/ESP during the flare (top) and detrended signals 
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with the previous method. The bottom panel shows the 
correlation coefficient as a function of time and period. 
Arrows are overlaid to show the phase difference between 
the channels. The cross-hatched area correspond to the 
so-called "cone of influence" , where results are influenced 
by edge effects. 

Not surprisingly, we find enhanced wavelet coherence 
during intervals la and lb for periods between 8 and 
32 s. This shows th at those fluctuations are broad- band, 
like for most QPPs (|Nakariakov and Melnikovll2009fl . We 
also find common oscillations during the impulsive phase 
in the 1-2 minutes range (lasting at least 2 periods). In 
both ranges of periods, the phase delay is close to — 7r/2; 
that is to say the EUV is leading the SXR channel by 
about a quarter of a period. During interval 2, we find 
common oscillations with periods of 3-5 minutes. This 
time, it is the EUV channel that lags the SXR channel 
up to 7r/2. This range of periods is reminiscent of the 
p- modes. 10 

We also find lagging and leading behaviors in wavelet 
coherence for other channels during intervals la and 

10 P-modes have bee n analyzed with the ESP SXR channel by 
Didkovsky ct al. (2011) for "quiet" (nonflaring) time periods. 



lb (not shown). In particular, the RHESSI 3-6, 6-12, 
12-25, and 25-50 keV channels show matching wavelet 
power in the 1-2 minute range with the ESP SXR chan- 
nel, and all lead this channel by about tt/2. As a cross- 
verification, the wavelet coherence analysis shows that 
they are nearly in phase with the ESP 18 nm channel in 
the 1-2 minute range. 

The Nobeyama 3, 9 and 17 GHz radio channels also 
show matching wave power with the ESP SXR and 18 nm 
channels in the 1-2 minute ranges; they are nearly in 
phase with ESP SXR, but lag ESP 18 nm by about -tt/2. 
For the 8-32 s-period range, their phase is different dur- 
ing intervals la and lb (« —tt/2 and ~ 0, respectively), 
but consistent with the analysis of the time lags pre- 
sented in Figured] 

To determine the spatial localization of the QPPs emis- 
sion, one must use images of the solar atmosphere. Un- 
fortunately, in most of the bandpasses (e.g. SDO/AIA, 
Hinode/XKF, RHESSI) the cadence of the image ac- 
quisition is too low. The only exception is the NoRH 
data at 17 GHz, which have 2 s cadence, sufficient to 
detect the sources of QPP. Figure 2] shows that a ra- 
dio source was present during interval la (upper left 
panel) and progressively moved to the West at 01:51. 
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It progressively faded, but was still visible during inter- 
val lb (lower left panel). Such moti on is not unusual, 
e.g. iZimovets and Struminskvl (|2009ft . A second source 
started to be resolved at 01:53:08; before, both sources 
overlapped within the instrumental resolution. From the 
comparison of the detrended light curves of both sources 
and that of the full Sun, we conclude that the fluctua- 
tions during both intervals la and lb are then emitted 
by distinct structures. 

4. DISCUSSION AND CONCLUSIONS 

We analyzed fluctuations on time scales « 8-32 s dur- 
ing the X2.2 flare of 15 February 2011. Their fine struc- 
ture was similar in EUV, radio, SXR and HXR channels 
during the impulsive phase, but with time delays up to 
~ 9 s. During the declining phase, the presence of fluc- 
tuations is less certain; they are only visible in the SXR 
channels, but still with time delays between channels. 

There have been many studies of time delays in fluctua- 
tions during flares ; for example between radio wavebands 
and HXR (e.g. Kaufman n et al.l 119831 : iCornell et al.l 
1984) or between HXR channels of different energy 
(Aschwandcn et al. 1997). To our knowledge, this is the 
first time that EUV fluctuations leading SXR fluctua- 
tions have bee n observed on these time scales. We note, 
however, that lEmslie et al.l (|1978l ) measured lags < 5 s 
in EUV radiation relative to HXR for pulses on similar 
time scales as those we analyze. 

Time shifts between clocks are of concern, but do not 
affect our general conclusions because there are signifi- 
cant time delays between channels of the same instru- 
ment, like the EUV and SXR channels of EVE/ESP. 

A key question is the emission mechanism at work in 
the different wavebands. There is no doubt that HXR are 
emitted through bremsstrahlung (thermal or most prob- 
ably non-thermal). But in the absence of observations 
with both high spectral and time resolution (« 1 s), we 
have no evidence of whether the SXR or EUV fluctua- 
tions occur in emission lines (through collisional excita- 
tion of the ions, mainly by thermal electrons) or in the 
continuum (implying bremsstrahlung). 

Three physical mechanisms can produce time de- 
lays in the range o f HXR energies (Aschwandcn 120041 : 
iHolmanet all 1201 1(1 : time-of-flight (TOF) dispersion of 
free-streaming electrons, magnetic trapping with the col- 
lisional precipitation of elect rons, and cooli ng of the ther- 
mal plasma (Neupert effect, Ncupcrt 19681). 

Electrons of higher energy have smaller TOF from loop 
top to the chromosphere, where they typically produce 
hard X-ray photons. Therefore, the TOF effect intro- 
duces 10-100 ms delays in the 25-50 keV light curv es as 
compared to 50-100 keV (Aschwandcn et al. 1995). As 
a consequence, TOF delays cannot explain our observa- 
tions: first, we observe larger delays than expected be- 
tween the HXR channels (« Is); second, the 50-100 keV 
channel should lead the 25-50 keV channel during inter- 
val lb. 

Due to magnetic mirroring between the loop foot- 
points, electrons of higher energy precipitate later be- 
cause they are les s likely to escape the magnetic trap. 
lAschwanden et all ((19971 ) finds typical delays of 1-10 s 
between HXR pulses at 200 keV and 50 keV and model 
them with trapping effects. The wis delay between the 
two HXR channels during interval lb is compatible with 



their results. The mismatch in the order of the channels 
in interval la could be explained by the effect of uncer- 
tainties, especially as this interval is shorter and contains 
fewer peaks for cross-correlation. 

Radio signals at 17 and 9 GHz are delayed with respect 
to the 25-50 keV signal by about 2 s during intervals la 
and lb. This is expected for trapped highly relativistic 
electrons (100-1000 keV) tha t produce g yro-synchrotron 
emission in those wavebands (|Dulklll98 5). The variation 
in delay from interval la to lb shown by HXR and radio 
channels coincides with different spatial locations of the 
radio sources (Figure |4|) . It is likely that the physical 
conditions (loop length, magnetic field) changed so that 
trapping times differ during those time intervals. As ex- 
plained above, results in interval la must be interpreted 
with caution, though. 

Cooling effects introduce delays in the peak times of 
channels more dominated by thermal effects (e.g. SXR) 
with respect to channels that pres ent a more impul- 
sive behavior (m icrowaves, HXR; e.g. [Aschwandcn] |2007l : 
Uain et ail 1201 If) . Cooling effects are compatible with 
short-period fluctuations in SXR channels lagging those 
in the impulsive HXR channels, as we observe. However, 
the same cooling effects are responsible for the flare peak 
to occur first in the GOES 0.5-4 A then in the GOES 1- 
8 A and in ESP 0-7 nm channels (see Figure[IJ . Paradox- 
ically, we observe exactly the reverse order in the short 
ti me scale s fluctuations (Figure [5]) . 

iDennisI (|1985l ) reports on nearly simultaneous (< 1 s) 
pulses in UV and HXR, as a result of chromospheric 
heating by precipitating electrons. We suppose that the 
EUV fluctuations during intervals la and lb also origi- 
nate from chromospheric heating, which explains the ob- 
served simultaneity with HXR channels (< 1 s during 
interval lb). 

Can the obser ved fluctuations be modulat ed by 
MHD waves? iVan Doorsselaere et al.l ([2011) used 
PROBA2/1YRA to observe oscillations with two dis- 
tinct periods during a flare: both « 8.5 s and « 75 s 
period oscillations appear in the Zr, Al and Lyman a 
channels. They interpreted the shorter periods in terms 
of the standing fast sausage mode and the longer pe- 
riods in terms of the standing slow sausage mode. It is 
tempting to interpret our observations similarly as stand- 
ing fast and slow sausage modes for the shorter (8-32 s) 
and longer (1-2 minute) periods respectively. However, 
apart from the similarity with the mode periods, our ob- 
servations cannot confirm or invalidate the modulation of 
the electron beam by MHD waves. Besides, MHD waves 
alone can hardly explain the observed delays between the 
signals in HXR, EUV and SXR channels. 

In summary, we suggest that the fluctuations observed 
during the impulsive phase (intervals la and lb) were 
produced by a beam of precipitating electrons, possibly 
modulated by MHD waves. When electrons reached the 
denser layers, they excited EUV as well as HXR emission 
in lower energy channels (12-25 and 20-50 keV). Higher 
energy electrons were delayed by trapping effects: the 
higher the energy, the longer the delay to escape the 
trap. Therefore, the fluctuations appeared a few seconds 
later in 50-100 keV HXR and then in 9 and 17 GHz radio 
channels. 

The case of the SXR channels is more puzzling. The 
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(full Sun) and NoRH detrendcd light curves for sources labeled "1" and "2" (right panel). 



delays in these channels compared to EUV and HXR 
channels could be due to the Neupert effect, although the 
relative delays between the SXR channels are in reverse 
order than what is expected in this case. This point 
certainly deserves more investigation and modeling. 
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time lags shown in Figure [2] 
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